Power-law statistics and stellar rotational velocities in the 

Pleiades 

J. C. Carvalho^l] R- Silva^'^Q J. D. Jr. do Nascimento^0 and J. R. De Medeirosfl 

^ Universidade Federal do Rio Grande do Norte, 
UFRN, Departamento de Fisica C. P. 16^1, 
Natal - RN, 59072-970, Brazil and 
'^Universidade do Estado do Rio Grande do Norte, 59610-210, Mossoro, RN, Brasil 

(Dated: March 4, 2009) 

Abstract 
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I. INTRODUCTION 



Some restrictions to the applicability of the statistical mechanics have motivated the 
investigation of the power-law or non-gaussian statistics, both from theoretical and experi- 
mental viewpoints. In this concern, the nonextensive statistical mechanics l| and extensive 
generalized power-law statistics [2'] are the most investigated frameworks. In the later one, 
recent efforts on the kinetic foundations of the /t-statistics leads to a power-law distribution 
function and the K-entropy which emerges in the context of the special relativity and in the 
so-called kinetic interaction principle (see Ref. [21]). Formally, the k- framework is based on 
K-exponential and ^-logarithm functions, which are defined by 



exp 



.(/) = (v/^+^+«:/)'/^ (1) 

ln«(/) = ^-j;^, (2) 



whereas the K-entropy associated with the ^-statistics is given by 

S^ = -J d'pf\n,f = -{\n,if)). (3) 

The expressios above reduces to the standard results in the limit k = 0. 

The Tsallis statistics has been investigated in a wide range of problems in physics ^. In 
the astrophysical domain, the first applications of this powerlaw statistics studied stellar 



polytropes [3|] and the peculiar velocity function of galaxy clusters {4]. More recently, Ka 



12|, 



niadakis statistics has also been studied in the theoretical and experimental context 
however the first application with a possiple connection with astrophysical system has been 
the simulation in relativistc plasmas. In this regard, the powerlaw energy distribuition pro- 
vides a strong argument in favour of the Kaniadakis statistics jsl. Therefore, in such systems 
where nonextensivity holds, the classical statistics may be generalized. In this sense, one of 
the most puzzling questions in stellar astrophysics in the past 50 years is that concerning 
the nature of the statistical law controlling the distribution of stellar rotational velocity, in 
spite of the large acceptance that stellar rotation axes have a random orientation [6;]. In 
the middle of the past century, Chandrasekhar and Munch were the first to derive ana- 
lytically the distribution of stellar projected rotational velocity, on the basis of a Gaussian 



^ For a complete and updated list of refences see http://tsallis.cat.cbpf.br/biblio.htm 
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distribution. In their approach, these authors first assumed a parametric form for a function 
f{v), where v is the true rotational velocity, then computed the corresponding distribution 
of the projected rotational velocity Vsini and finally adjusted a set of stellar parameters 
to reproduce the Vsini measurements. Two decades later, Deutsch claimed that the 
distribution of stellar rotational velocities should have the form of a Maxwellian-Boltzmann 
law. Nevertheless, a number of studies have shown a clear discrepancy between theory and 
observations, where observed distributions are not fitted by a Gaussian or Maxwellian func- 
tion with a good level of significance. A Gaussian or Maxwellian distribution that fits the 
fast rotators fails to account for low rotation rates. On the other hand, a fit to slow rotators 
fails to explain the rapidly rotating stars j^, [lo| . 

Rotation is one of the most important observable in stellar astrophysics, driving strongly 
the evolution of stars, providing also valuable informations on stellar magnetism, mixing 
of chemical abundances in stellar interior, tidal interaction in close binaries, and engulfing 
of brown dwarfs and planets. In addition, if the present value of the rotational velocity 
of stars at a given evolutionary stage reflects the original angular momentum with which 
they were formed, the behaviour of the distribution of rotational velocity may also be used 
to study some of the characteristics of the physical processes controlling star formation. 
Early studies on the nature of the statistics controlling the distribution of stellar rotational 
velocity were based on Vsini measurements with poor precision, which, admittedly, can 
lead to systematic errors on the final analyses for low Vsini values. Here, we show that the 
question of the nature of the distribution of stellar rotational velocity, at least for low-mass 
stars in the Pleiades open cluster, is not simply a question of which mathematical function 
model is used, but it depends primarily on the statistical mechanics applied, which should 
be general enough to take into account the changes in rotation with time. 

In this work we have investigated the effects of the powerlaw statistics on the observed 
distribution of projected rotational velocity measurements of stars in the Pleiades open 
cluster by considering a ^-distribution functions. In this regard, it is worth emphasis that 
earlier study based on the Tsallis statistics has been developed in Ref . [l^ . However, in the 

present work, to study the effects of the powerlaw statistics on the observed distribution in 

i 

the Pleiades open cluster we use the more recent generalization of Kaniadakis [J] and, for 
completeness, we compare the results with the ones obtained in the context of the Tsallis 
statistics. This paper is organized as follows. In Sec II, based on the basic formalism 
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presented in Ref. js], we present a generalization of the rotational velocity distribution in 
the spirit of Kaniadakis statistics. A brief discussion on the stellar sample is made in Sec. 
III. Our main results are discussed in Sec. IV and we summarize the main conclusions in 
Sec. V. 



II. K-DISTRIBUTION FUNCTION 

As well know, a large portion of the experimental evidence, as well as some theoreti- 
cal considerations supporting Tsallis and Kaniadakis proposal involves a non-Maxwellian 
(powerlaw) distribution function associated with the thermostatistical description from the 



variety of the physical systems [ll|, ll2|. In Tsallis framework, the equilibrium velocity q- 
distribution may be derived from at least three different methods, namely: (i) through a 
simple nonextensive generalization of the Maxwell ansatz, which is based on the isotropy of 



the velocity space 13[; (ii) within the nonextensive canonical ensemble, that is, maximizing 



Tsallis entropy under the constraints imposed by normalization and the energy mean value 
and (iii) using a more rigorous treatment based on the nonextensive formulation of the 



Boltzmann if-theorem [15|. Here, we revisit the first method by considering the Kaniadakis 
statistical which is based on exp^ and ln«; given by Eq. ([1]), where / is a function of random 
variables that includes the standard exponential as the limiting case when k ^ 0. 

It is widely known that Deutsch [8'] has considered the distribution function for the 
magnitude of a vector that has random orientation. For this, it is required to find the 
distribution function of a positive scalar u;, which is the magnitude of a vector Q. We 
assume that the distribution of u is isotropic. We also assume that if it is decomposed into 
components along Cartesian axes, the distribution of any component is independent of the 
other components. Deutsch has defined f2 as the non-dimensional quantity uoj, where j is a 
parameter with the dimension uj~'^ , so 

= QJ+ Vty] + VtS (4) 

The probability that VL.j. lies in the interval VL^ + (ifi^,.], Vty in [VLy \ Vty + dVly] and VL^ in 
Viz + dVlz] is then 

Fisi)^^! = f{n,)f{ay)f{n,)dn,dnydn, (5) 



4 



where fl = ^yTf^^T-TTT^'+'ni and F{Q) is the standard Maxwelhan distribution function 



F{n) = ^n^exp{-n^) (6) 

Jit 



Let us now consider the arguments given in Refs. [l3|,ll6(. One can modify the basic 
hypothesis of statistical independence between the distributions associated with the com- 

- n n 

ponents of Q, based on the ^-statistics. As pointed out in and pj[, the independence 
between the three velocity components does not hold in systems with long-range interaction, 
or statiscally correlated, where the power-law statistics character is observed. Taking such 
arguments into account, the generalization for Eq. (6) in the light from Kaniadakis the 
framework reads, 

F{n)d^n = exp,(ln, f{n,) + In, fi^y) + In, /(a))dfi,c/fi,cin„ (7) 

where the K-exponential and /t-logarithm are given by identities (1) and (2). In particular, 
in the limit k = the standard expression (JSj) is recovered. Note also that lnK(exp,(/)) = 
exp,(ln,(/)) = /, and ^^^^ = — ^ are satisfied. Therefore, the partial differentiation 
of the K-ln of with respect to Qi leads to 

dln^F d 



or, equivalent ly. 



■(In, f{n,) + In, /(^],) + In, /(^],)), (8) 



F'^ + F-'^ IdF _ 1 d 
2F xd^ ~ a rfa 

where x = ^/^i + + ^'t- Now, by defining 

, , F^ + F-" IdF 

we may rewrite ([9]) as 

The second member of the above equation depends only on with i = x,y,z. Hence, 
equation ( ITTi) can be satisfied only if all its members are equal to one and the same constant, 
not depending on any of the velocity components. Thus, we can make = ~2/o"^, where 
the parameter cr, is the width of the /t-distribution, leading to 

1 (9 2 

-(ln,/,) = --. (12) 



Hence, the solutions of Eq. (fT2|) for f{Qi) is given by the ^-distribution 

f{Q,) = exp^{-Q^/al) (13) 

From f[T^ we see that the Gaussian probabihty curve is replaced by the charactheristic 
powerlaw behavior of Kaniadakis framework and, as expected, the limit k, = recovers the 
exponential result. Note also that for any values of the K-parameter, the powerlaw (fT3l) does 
not exhibits a cut-off in the maximal allowed rotational velocities, and it is straightforward 
to show that F{Q) is given by 

Fm = exp^i-n'/al) (14) 

The probability of finding uj in the interval [Q, Q+dQ] is determined by ^(f^) = / f{Q)d^fl 
which leads to a power law that belongs to the same class of power law as given in Eq. f|T3l) . 
i.e., 

^{n)=Ann'exp,i-n'/al). (15) 

Here, it is worth mentioning that the standard distribution of the true rotational velocity V 
for a star sample is F{V) ~ V'^exp{—V^). As shown by Deutsch j^, the standard observed 
distribution of the projected rotational velocity Vsini, for a random orientation of axes, 
must be given by (f){y) ~ y exp{—y'^) [l^], with y = Vsini. Henceforth, the ^-distribution 
(pniy) should reproduce the standard one, in the same way as Fi^{v) recovers F{V) in the 
K = limiting case. Therefore, by considering this arguments, we introduce the following 
distribution function for the observed stellar rotational velocities 

= yexp^i-y'^/al). (16) 



III. THE STELLAR SAMPLE 



The rotational velocities V sin i used in the present analysis were taken from the rotational 



survey for the Pleiades stars carried out by Queloz et al. [18|. We have selected 219 stars 
from the original sample given by those authors. All the selected objects are low-mass 
stars and provide a complete and unbiased rotation data set for stars in the B ~ V range 
(0.4—1.4) corresponding to an effective temperature range from 4000 K to 6000 K and a mass 
range from 0.6Mq to 1.2Mq. For a complete discussion on the observational procedure. 



TABLE I: Best values of the parameters of Kaniadakis (k and a^) and Tsallis {q and aq) distribution 
determined using Kolmogorov-Smirnov test for the rotational velocity of stars in the Pleiades 
cluster. 



A{B - y) N K C7k Pmax Q CFq Pfnax 



0.40-1.40 217 0.4461^:^^^ 7.81 0.28 1.334+^;^^° 6.93 0.23 



calibration and error analysis the reader is referred to Queloz et al. [18| . We should observe 
that individual errors in V smi measurements are better than about 1 km/s and should not 
play a significant role on the observed distributions. 

If we plot the observed distribution we note that there are two stars with velocity (105 
km/s and 160 km/s) too far from the peak of the distribution (5-8 km/s). We have, 
therefore, excluded these stars with exceedingly high value of V smi from the sample. 

In order to avoid biases due to arbitrary choices of bin range when constructing the 
frequency histograms, we have decided to study the observed cumulative distribution of the 
rotational velocities, V sini, and compare it with the integral of the probability distribution 
function in (fT6|) . The normalized cumulative distribution is given by 

/o°°2/exp^(-2/V(T2)d|/' 

IV. RESULTS 

Parallel with the Kaniadakis distribution we have also fitted the data to the Tsallis 
distribution in order to evaluate which one would best fit the observations. The Tsallis 
generalized maxwellian is a two parameter [q and cTq) function given by 

^21 1/(1-9) 
72 



'1 



{y) = y 



(18) 



In the limit q = 1 the standard Maxwellian is recovered. 
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FIG. 1: Observed distribution (histogram) of the rotational velocity of the stars in the Pleiades 
open cluster. The curves represent the best fitted maxwellian (dashed line), Tsallis (thin line) and 
Kaniadakis (thick line) distribution. The fitting parameters are in Table 1. 

To calculate the best values of the distributions parameters for the observed cumulative 
distribution we used the Kolmogorov-Smirnov statistical test. 

The distribution functions were used to fit the observational data, to obtain the best 
(pK^y) and 4>q{y), giving the best k and g- value together with the best 0"^ and cTg for the 
corresponding distribution. The results are shown in Table 1. We can clearly see that the 
Vsini distribution of the Pleiades stars do not obey a standard Maxwellian function since 
the values of k and q are significantly different from and 1, respectively. 

Fig. [Ushows the best fits for the histogram of the observed distribution of V sin i according 
the results in Table 1. The k and q-Maxwellian functions are represented, respectively by the 
thick [k, = 0.45) and thin {q = 1.33) lines. The dashed line represents standard Maxwellian 
function. The distribution of observed V sin i is without a doubt more adequately fitted 
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FIG. 2: As in Fig. [T] but for the logarithm of the distribution function divided by Vsini as a 
function of (Vsini)'^. 

by either the k or the q-MaxwelUan function. This is more noticeable in Fig. [2] where we 
have plotted the logarithm of the distribution divided by Vsini, that is log{(j){y) /y), as a 
function of (\^sinz)^ so that the standard Maxwellian is represented by a straight (dashed) 
line. The Kaniadakis distribution is represented by the thick line while Tsallis distribution 
by the thin line. We observe that none of the distribution fit well the high velocity end 
of the observed data. Although the Kaniadakis function fits the data slightly better than 
Tsallis function, the difference may be regarded as marginal as indicated by the values of 
the maximum probability of the Kolmogorov-Smirnov statistical test (Table 1). 

Finally, in Fig. [3] we present the behaviour the parameter a, representing the width of 
the two non-Maxwellian, as a function of the parameter k and q. It is clear that, at least 
for the present stellar sample, the standard Maxwellian (k = or g = 1) is in the rejection 
region, i.e., outside the curve which delineates 0.05 significance level. 
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FIG. 3: Rejection region (outside curves) of the null hypothesis that the V sin i distribution is drawn 
from the ^-distribution (upper panel) at 0.05 and 0.15 confidence level and from g-distribution 
(lower panel) at 0.05 and 0.10 confidence level. The dot (•) represents the maximum probability 
for the pair (k - a^) and {q - (jq). 

V. CONCLUSIONS 

In this work we have used non-gaussian statistics to investigate the observed distribution 
of projected rotational velocity of stars in the Pleiades open cluster. We have studied 
in details the Kaniadakis distribution and shown that it fits more closely the observed 
distribution than the standard maxwellian. A comparison with the Tsallis statistics shows 
that, at least when the observed distribution presents an extended tail, as it is the case of 
rotational velocity of stars in the Pleiades, both distributions give equivalent, though not 
entirely satisfactory results. 

As discussed in Sec. IV, the best fits for the histogram of observed distribution are non- 
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gaussian with k = 0.45 and q = 1.33 for Kaniadakis and Tsallis parameter, respectivily. 
In particular, we emphasize that the resuh of the g-parameter is consistent with the upper 
hmit q < 2 obtained from several independents investigations in the quantum limit [l^ and 



in the non-quantum limit 



20|. 



Finally, it is worth mentioning that the the best fits in k, and q can be recalculated by 
considering a more robust stellar sample. In this respect, the stellar radial velocity of a 
sample of open clusters are being studied. This issue will be addressed in a forthcoming 
communication. 
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